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We theoretically investigate the physics of interacting Bose-Einstein condensates at equilibrium in a weak 
(possibly random) potential. We develop a perturbation approach to derive the condensate wave function for 
an amplitude of the potential smaller than the chemical potential of the condensate and for an arbitrary spatial 
variation scale of the potential. Applying this theory to disordered potentials, we find in particular that, if 
the healing length is smaller than the correlation length of the disorder, the condensate assumes a delocalized 
Thomas-Fermi profile. In the opposite situation where the correlation length is smaller than the healing length, 
we show that the random potential can be significantly smoothed and, in the mean-field regime, the condensate 
wave function can remain delocalized, even for very small correlation lengths of the disorder. 
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I. INTRODUCTION 

Ultracold atomic gases are currently attracting a lot of atten- 
tion from both experimental and theoretical viewpoints. Tak- 
ing advantage of the recent progress in cooling and trapping 
of neutral atoms [1], dilute atomic Bose-Einstein condensates 
(BECs) 1 2] and degenerate Fermi gases (DFGs) S H] 
are now routinely produced at the laboratory. Using vari- 
ous techniques, space-dependent potentials can be designed 
almost on demand in these systems. For example, one can 
produce periodic p,|81, quasi-periodic [9, 10, 1 IJ, or random 
potentials 

iilHlllQlMIll by using optical means. For 
these reasons and due to unique control and analysis possibil- 
ities, ultracold gases constitute a favorite playground for re- 
visiting standard problems of condensed matter physics (CM) 

Most current experiments with BECs lie in the mean-field 
regime where the Bose gas is described by a single wave func- 
tion, ip, governed by the (nonlinear) Gross-Pitaevskii equation 
ll22ll . Due to the interplay between the kinetic energy term and 
the interaction term, it is usually difficult to derive the exact 
solution of this equation. The importance of interactions can 
be characterized by the healing length, ^, which defines the 
typical distance below which the spatial variations of sig- 
nificantly contribute to the energy of the BEC, via the kinetic 
energy term ll22ll . In the Thomas-Fermi regime (TF), i.e. when 
^ is significantly smaller than the typical variation scale, CTr, 
of the potential, V{r), to which the BEC is subjected, the ki- 
netic term is negligible and the BEC density simply follows 
the spatial variations of the potential': 



(1) 
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' This is standard in the case of a harmonic confinement, V{r) = 
muP'r'^ jl. Although, there is no intrinsic typical variation scale, one can 
define ctr as = fi, i.e. iTr = Ltf, the usual TF half size of 

the condensate and the validity of the TF regime reads ^ <g Ltf- For pe- 
riodic, quasi-periodic or random potentials, ctr is the spatial period or the 
correlation length (see section|lIl]for details). 



In the opposite situation (^ > ctr), the kinetic term should be 
taken into account and the exact BEC wave function usually 
cannot be found analytically. 

Besides a general interest, the question of determining the 
BEC wave function for an arbitrary ratio CTr/^ has direct ap- 
plications to the case where V{r) is a random potential. The 
physics of quanturn systems in the presence of disorder is 
central in CM | 23l |24 l25ll . owing to unavoidable defects in 
'real-life systems' . One of the major paradigms of disordered 
quantum systems is due to Anderson who has shown that the 
eigenstates of single quantum particles in arbitrary weak ran- 
dom potentials can be localized, i.e. ip shows an exponential 
decay at large distances^ i26ll . Recent experiments have stud- 
ied the onset of strong or weak localization effects of light 
waves ll27i l28ll and microwaves ll29l [soll . Ultracold matter 
waves are also widely considered as promising candidates to 
investigate Anderson localization in random |31, 32, 33] or 
quasi-random structures 13 ll l34tl and more generally to 
investigate the effects of disorder in various quantum systems 
(for a recent review, see Ref. ll35ll and references therein). It is 
expected that the dramatic versatility of ultracold gases would 
allow us for a direct comparison with theoretical studies of 
quantum disordered systems. 

A key peculiarity of BECs is that interactions usually can- 
not be neglected and interaction-induced delocalization can 
compete with disorder-induced localization effects |15, 1^ 
[TtIi . Generally, the interplay between the kinetic energy, the 
interactions and the disorder is still a o pen question that has 
motivated many works |36, 37, 38, [39ll40ll . It is clear from 
Eq. ([T]i that, in the TF regime (ctr ^ 0> where the interaction 
forces the wave function to adapt to the random potential, a 
BEC will not localize. Indeed, if V{r) is a homogeneous ran- 
dom function^ |41], so is the BEC wave function, tp, which 
therefore, cannot decay at large distances. This has been con- 



In ID and 2D systems, all eigenstates are usually localized while in 3D, 
they are localized below the so-called mobility edge. 
^ In this context, the term 'homogeneous' means that all local statistical 
properties of the random potential are independent of the position. 
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firmed in recent experiments ifTsi [161 [TtIi . The question thus 
arises as to understand whether, as a naive transcription of 
the loffe-Regel criterion ll42ll would suggest, localization can 
happen when ctr < ^. 

In this paper, we show that this criterion is actually not 
sufficient for BECs at equilibrium if the interactions are non- 
negligible (i.e. if ^ <C L, where L is the size of the system). 
We indeed show that interaction-induced delocalization still 
overcomes localization effects even when ^ ^ CTr. In fact, 
due to the smoothing of the random potential f^^, the effect 
of disorder turns out to be reduced when ^/ctr increases. 

In the following, we develop a general formalism based 
on perturbation theory (see section to determine the BEC 
wave function in any given weak potential, V{y), for an arbi- 
trary ratio ctr/^. We find that the BEC density, j^/ip, is still 
given by Eq. ([TJ, except that the potential V^(r) has to be re- 
placed by a smoothed potential, V{t). We derive an exact 
formula for the smoothed potential up to first order in the per- 
turbation series. We then apply our results to the case where 
V{y) is a ID homogeneous random potential (see section HiHi 
and derive the statistical properties of the smoothed random 
potential, V{y). From this, we conclude that an interacting 
BEC remains delocalized, even for ^ ^ ctr (if ^ ^ L). 



II. SMOOTHING EFFECT IN INTERACTING 
BOSE-EINSTEIN CONDENSATES 

Consider a low-temperature Bose gas in d dimensions 
with contact atom-atom interactions, gdY)5^''-\v), where gdx) 
is the d-dimensional interaction parameter. In 3D geome- 
tries, (73D = 47rfi,^asc/TO, where is the scattering length 
ll22ll . and m is the atomic mass. Low dimensional ge- 
ometries (ID or 2D) can be realized in ultracold atomic 
samples using a tight radial confinement, so that the radial 
wave function is frozen to zero-point oscillations in the form 
where rj^ is the radial coordinate vector. In this 
case, gd-D — Usd J drj^ |i^*[(rj^)|'*. For instance, one finds 
giu = 2hLU±asc, for a 2D harmonic radial confinement of 
frequency uj±. In addition, the Bose gas is assumed to be 
subjected to a given potential, V{r), with a typical amplitude 
K and a typical variation scale ctr. Possibly, the potential, 
V{r), may have various length scales. In this case, we as- 
sume that (Tr is the smallest. Assuming weak interactions, i.e. 
—2/d-i ^ mgdo/fi^, where n the mean density |44, 45], we 
treat the BEC in the mean-field approach 1221] and we use the 
Gross-Pitaevskii equation (GPE): 



2m 



Vir)+gdDmr)\' 



V^r) , (2) 



where /i is the BEC chemical potential, and where the wave 
function, ip, is normalized to the total number of condensed 
atoms, / dr |-0(r)p = A^. Note that ip minimizes the A^-body 
energy functional so that ^jj is necessary a real function (up 
to a non-physical uniform phase). In ID and 2D geometries 
and in the absence of trapping, no true BEC can exist due to 
significant long- wavelength phase fluctuations [4^. In this 



case, no macroscopic wave function, tp, can be defined. How- 
ever, because density fluctuations are strongly suppressed in 
the presence of interactions, the Bose gas forms a quasicon- 
densate ['46^ and the density, n, can be treated as a classical 
field. It turns out that ^/n is governed by Eq. (|2]i. Therefore, 
even though we only refer to BEC wave functions in the fol- 
lowing, our formalism also applies to quasicondensates, after 
replacing by ^/n. 



A. The Thomas-Fermi regime 

In the simplest situation, the healing length of the BEC is 
much smaller than the typical length scale of the potential 
(C ^ ctr)- Therefore, the kinetic energy term in the GPE (|2]i 
is small and the BEC density, simply follows the spatial 
modulations of the potential: 



|^(r)|2 = [^-y(r)]/5,D ^ox^i>V{v) 



and \il){v)\'^ = 



otherwise. 



(3) 



This corresponds to the TF regime. Note that for 14 <C /i, one 
has 



2fi 



(4) 



with ipo = ^/Ji/gdD being the BEC wave function at V{y) = 
0. Therefore, the BEC wave function itself follows the modu- 
lations of the potential V{y). 



B. Beyond the Thomas-Fermi regime: the smoothing effect 

The situation changes when the healing length is of the or- 
der of, or larger than, the typical length scale of the potential 
(^ > ctr)- Indeed, the kinetic contribution limits the small- 
est variation length of the spatial modulations of a BEC wave 
function to a finite value of the order of the healing length 
||22|] . Therefore, the BEC can only follow modulations of the 
potential on a length scale typically larger than ^ and Eq. (O 
no longer holds. 

For a weak amplitude of the potential"^, we can use perturba- 
tion theory techniques. We thus write the BEC wave function 
as -0(1") = "00 + 5'4){y) where we assume that Sip <C 'ipo, and 
V^o is the zeroth-order solution of the GPE (|2]i: 



^■00 



2m 



(5) 



Since the BEC is homogeneous at zeroth-order, one has 0o = 
yijJgdD- Then, the first order term of the perturbation series 
is given by 



2m 



V\5^P) - [/X - 3gdD0o] H = -"^'(r)Vo 



(6) 



^ A precise condition for the validity of the perturbative approach will be 
given later [see Eq. )16H . 
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Since ^ ~ Sgdoipo = ~2/x, we are left with the equation 



2ix 



(7) 



where ^ = h/ ^/2m/i is the heahng length of the BEC. We 
straightforwardly find the solution of Eq. (|7]l, which reads 



j dr' G(r - r' 



2fi 



(8) 



where G(r) is the Green function of Eq. O, defined as the 
solution of 



f2 



G{r) (r) 



or equivalently, in Fourier space 



1 



G(k) = l/(27r)'*/2 , 



(9) 



(10) 



where, G(k) — ^^^^^^^ / dr G(r)e~''' '" is the Fourier trans- 
form of G. In contrast to the case of single particles, the Green 
function, G(k), has no singularity point so that the perturba- 
tive approach can be safely applied for any wave vector k. 

The explicit formula for the Green function, G, depends on 
the dimension of the system. After some simple algebra, we 
find 



in ID, 



G(z) 



in 2D, G{p) = 



in 3D, G(r) = 



J_cx ( ^ 
1 



<2 



rKn 



1 



2^2 |r I 



cxp 



aV2 



(11) 

(12) 
(13) 



where Kq is the modified Bessel function. Finally, up to first 
order in the perturbation series, the BEC wave function reads 



2ii 



with 



V{r) ^ / dr'G(r')V^(r-r') 



(14) 



(15) 



Interestingly enough, the Green function in any dimension 
shows a exponential decay, with a typical attenuation length, 
^, and is normalized to unity^, J dr G(r) = 1. Therefore, 
G(r) can be seen as a smoothing function with a typical width 
^. Indeed, it should be noted that Eq. (fl4] i is similar to Eq. (|4|i, 
except that the potential V{r) which is relevant in the case 



^ This property follows directly from the definition jlO) of the Green func- 
tion. Indeed, / dr G(r) = (27r)<*/2G(k = 0) = 1. 



^ <C (Tr, changes to the potential V{r) for ^ > ctr. The poten- 
tial V{r) is a convolution of V{r) with a function which has a 
typical width ^ and thus corresponds to a smoothed potential 
with an amplitude smaller than K. In addition, if ctr corre- 
sponds to the width of the correlation function of a random 
potential, V, the correlation length of the smoothed random 
potential, V, is of the order of max(crR,^) [for details, see sec- 
tionHni''. 

Note that, for^ <C ctr, G(r) can be approximated by (5^'^^(r) 
in Eq. ( fTSl l. and V{y) ~ V{y). We thus recover the results of 
section UTaI valid for the TF regime. 

The validity condition of the perturbation approach directly 
follows from Eq. ( fT4l l: 



V{v) < 11 



(16) 



Note that if ^ ^ ctr, the potential can be significantly 
smoothed so that the above condition can be less restrictive 
than the a priori condition, V{v) <C ^i- 

The results of this section show that the potential, V{y), can 
be significantly smoothed in interacting BECs. We stress that 
this applies to any kind of potentials provided that ^ <C i and 
V{y) <C /i. In the next section, we present an illustration of 
the smoothing effect in the case of a random potential. 



III. APPLICATION TO A TRAPPED INTERACTING 
BOSE-EINSTEIN CONDENSATE IN A ID RANDOM 
POTENTIAL 

A. Trapped ID Bose-Einstein condensate in a random 
potential 

In this section, we consider a ID Bose gas subjected to a 
weak homogeneous random potential, V{z), with a vanishing 
average value ({V) = 0), a standard deviation, T4, and a spa- 
tial correlation length, ctr, significantly smaller than the size 
of the system. In addition, we assume that the gas is trapped 
in a confining harmonic trap^, V[,(z) = muJ^z^ 12 as in almost 
all current experiments on disordered BECs iflil ITsl [161 [TtIi . 
We consider a situation such that hio ^ ngi^ <C H^n'^/m, 
i.e. the Bose gas lies in the mean-field regime, and in the ab- 
sence of disorder, the interactions dominate over the kinetic 
energy^. The situation mimics the experimental conditions of 



* In contrast, for example in the case of a deterministic periodic potential, 
V{z) = VrCos (kz), the vaiiation scale, ctr = 27r/fc, corresponds to the 
period of the potential, and we find V{z) = ^-j^^i^^ ■ The potential is 

indeed smoothed as the amplitude of V is smaller than that of V. Never- 
theless, the period of the smoothed potential, V, is the same as that of the 
bare potential, V. 

^ All results also apply if there is no trapping. In this case, all zeroth-order 
terms simply do not depend on z. 

* This corresponds to the usual TF regime for confined BECs in the absence 
of disorder | 22]. However, no restriction is imposed for the ratio ctr/^, so 
that the BEC can be out of the TF regime as defined in section|l] 
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Figure 1 : (color online) Example of the realization of a speckle ran- 
dom potential with ctr ~ 10~^Ltf. 



The presence of the harmonic confinement in- 
troduces a low-momentum cut-off for the phase fluctuations 
so that the ID Bose gas forms a true condensate at low tem- 
peratures ll44ll45ll . In this case, the BEC wave function is 



"00 = y/Tdz)Jg 



ID 



(17) 



where ^JLQ{z) = /i — is the local chemical potential. 

This corresponds to an inverted parabolic density profile with 
a half-size, L^f — y^/i/mtJ^, where the chemical potential 

is ^ = = (M^i9i^^^'^^ with ; = y?^7^ being the 

extension of the ground state of the harmonic oscillator 

As Ltf ^ ('^, (Tr), it is legitimate to use the local density ap- 
proximation (LDA) [22], i.e. in a region significantly smaller 
than Ltf, the quantities ifj^ and can be considered as uni- 
form. We can thus directly apply the results of section HTbI 
From Eqs. (fT4ll-(fT7]i, we immediately find that 



n{z) ~ no(z) 



V{z)_ 

giD 



where 



V{z) 



dz' 



exp 



V2Uz) 



V{z - z') 



(18) 



(19) 



is the smoothed potential, with £,o{z) — ft/ \/2m^Q{z) being 
the local healing length. The density profile of the BEC is thus 
expected to follow the modulations of a smoothed random po- 
tential. 

Note that the total number of condensed atoms is = 
/ dz \^/nQ{z) + S^p\'^ ~ / dz (710(2) - ^^(2)/5id) up to first 

order in V / iJ.. Since {V) — 0, one has (A^) ~ Jdz no{z), 
owing to the assumed self-averaging property of the potential 
ll4lll . In addition, we have fi = fij^. 

We now compare our predictions to the exact solutions of 
the GPE ^ as obtained numerically. For the sake of concrete- 
ness, we consider a speckle random potenti al 14711 s imilar to 
the one used in the recent experiments ifTil Il5l Il6l [TtIi (see 
Fig.[T]i. Briefly, a speckle pattern consists in a random inten- 
sity distribution and is characterized by its statistical proper- 
ties. First, the single-point amplitude distribution is a negative 
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Figure 2: (color online) Density profiles of a BEC confined in 
a combined harmonic plus random potential (Vr = O.ljj,, ctr = 
7.5 X 10^'^I/tf). The solid (red online) line corresponds to the nu- 
merically computed BEC wave function; the dashed (green online) 
line is the TF profile in the absence of disorder; and the black dotted 
line is a plot of the disordered TF profile [Eq. a) Case where 
the healing length at the trap center, ^, is smaller than the correla- 
tion length of the random potential: ctr/^ — 10. In this case, the 
density profile follows the modulations of the random potential ac- 
cording to Eq. ((BJ. b) Opposite situation: (Tr/^ ~ 0.5. In this case, 
the BEC density profile, obtained numerically, significantly differs 
from Eq. but can hardly be distinguished from Eq. (Tsj [also 
plotted in Fig.[2j)) as a dotted (purple online) line]. The inset shows 
a magnification of the plot in a very small region of the BEC. 



exponential 

r[v{z)] 



exph(y(z) + K)/K] 
IKI 



(20) 



and = otherwise , 

corresponding to the average value {V) = and the standard 
deviation AV = x/{[V{z) ~ (y)]2) = |K|. Second, the spa- 
tial correlations are characterized by the autocorrelation func- 
tion C(Az) = (I^(Az)V"(0)) which correlation length is de- 
noted Ctr and can be chosen at will lfl7ll47l] . For the numerical 
calculations, we numerically generate a ID speckle pattern us- 
ing a method similar to the one described in Ref. [48] in ID 
and corresponding to the correlation function 



C(A2) = sinc(\/3Az/\/2CTR) 



(21) 



where sinc(a;) — sin(a;)/x. For the sake of simplicity, it is 
useful to approximate C(z) to a Gaussian function (see for 



5 



example section llllBl l. Up to second order in ^zja^, we have 
C(Az)c.K'exp(-AzV2a2). 

Numerical solutions of the GPE (|2]i are presented in Fig. |2] 
for two values of the ratio CTr/^, where ^ is the BEC healing 
length at the trap center. In the first case (Fig. |2^), we have 
^ ^ (Tr, and the density simply follows the modulations of 
the bare random potential, according to Eq. (O. In the second 
case (Fig. I2J5), we have ^ > ctr, and as expected, the BEC 
wave function does not follow the modulations of the bare 
random potential V{z) but actually follows smoother mod- 
ulations of the smoothed potential V{z). Figure |2]3 (and the 
inset) shows that the numerically computed density can hardly 
be distinguished from Eq. ( fTsT l. This supports the validity of 
our approach. 

B. Statistical properties of the smoothed random potential 

It is useful to compute the statistical properties of the 
smoothed random potential V(z) as they will be imprinted on 
the BEC density profile according to Eq. ( fTsT l. From Eq. ( fTsT i, 
we immediately find that, (i) V(z) is a random homogeneous 



potential, (ii) the average value of V vanishes, 

{V) = (V^) = , (22) 
and (iii) the correlation function of V is given by 

a(Az)- y" AuAvC\^z + (v-u)\G,{u)G,+i^,{v) , (23) 

where C(Az) = (V"(Az)F(0)) is the correlation function of 
the bare potential V and Gz (u) is given by Eq. ( fTTT i with ^ re- 
placed by ^q{z). In the following, we assume that Az ^ Lj^ 
so that Gz — Gz+Az and we omit the subscripts. Assuming 
for simplicity a Gaussian correlation function for the bare ran- 
dom potential, C(Az) ~ exp(— Az^/2f7^), we find after 
some algebra 

c(Az) = v;^i](|,^) ' ^'^^ 

with 



S(ctr, Az) 



-a^(^l- 2al - V2 Az) exp (wl + V2 Az) erfc 



2aR 



(25) 



-^CTr (^l-2al + V2 Az) exp (wl - V2 Az) erfc 



^\ . / - \/2 Az 



2ctr 





Figure 3: (color online) Left: Plot of the correlation function 
E ^1^, ^). Right: Width at l/^/e of the normalized correlation 

function S 



where CTr = ctr/^q, Az = Az/^o and 
erfc {x) ~ dt e~* is the complementary er- 

ror function. The correlation function S ^) is plotted 
in Fig. [3] 



This function S clearly decreases with CTr/^o^ indicating the 
onset of an increasing smoothing effect. At Az = 0, we have 
a simple asymptotic expression for CTr 3> ^o^ 



E(aR/eo,0)^l 



CTr > Co 



(26) 



So, as expected, I]((Tr/CoiO) — > 1 as CTr/^o — ^ oo. i-e- the 
random potential is hardly smoothed. For CTr ^ ^o- 



S(fTR/eo,0) 



%/7r ctr 
2 ^0' 



o-R < ^0 



(27) 



So, J^{(Tfi/^o,0) ^ as (Tr/^o 0, i.e. the amplitude of 
the smoothed random potential is significantly reduced com- 
pared to the amplitude of the bare random potential. Gen- 
erally speaking, from Eq. ( l24b . we have {V'^) — C(0) — 
V^Y.{aR/£,o,0). It follows that {V"^) is an increasing func- 
tion of (Tr/^o and that {V'^) < V^. This is consistent with the 
idea of a smoothing of the random potential. 
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In addition, the correlation length, ctr, of the smoothed ran- 
dom potential, V, is given by the width at of the func- 
tion Az S(ctr/^o, A-^/fo)- At (Tr the smoothing 
is weak and CTr ~ ctr. For ctr < i^o, the smoothing is signif- 
icant, so that Ctr saturates at CTr ~ as expected. Roughly 
speaking, we have CTr ~ ma.x{a^, ^q) [see Fig-E). 



C. Effect of disorder in interacting Bose-Einstein condensates 

We finally discuss the properties of the BEC wave func- 
tion in the presence of disorder It follows from Eq. ( fTSI l that 
the BEC density follows the modulations of a random poten- 
tial V. In the TF regime <C ctr), V :^ V, while when 
^ > (Tr, y is smoothed. Since y is a homogeneous random 
potential, there is no decay of the wave function. In particu- 
lar, Anderson localization does not occur, even for ^ ::t> ctr. 
In the case when ^ > ctr, it turns out that the BEC density is 
actually less affected by the random potential than in the TF 
regime {(_ <C ctr). This is in striking contrast with the case of 
non-interacting particles where localization effects are usually 
stronger at low energy 1I4TI1 . 

More quantitatively, using the statistical properties of the 
smoothed random potential, V, one can easily compute the 
fluctuations An{z) = ^ — 710(2;)]^) of the BEC den- 
sity around the average value riQiz) = \\i — muP' l'2\l g^. 
From Eq. we find An^ ~ C{^)lg^^. Note that At? de- 
pends on the displacement from the trap center through the 
dependence of ^0 on z. At the trap center, we find 



An, = ^\/S(VP)- 

5lD 



(28) 



We recafl that ^ = ^o(O) = h/^/hnyL is the BEC heaHng 
length in the trap center 

We have numerically extracted the fluctuations of the den- 
sity in the trap center, according to the formula An^ ~ 

\J -rtn Hz) - no{z)]\ This provides a good 

estimate of Auc as (,o{z) changes by less than 3% in the range 
[— Ltf/4, +Ljp/4:]. As shown in Fig.|4] the numerical results 
perfectly agree with Eq. d28T l over a large range of the ra- 
tio ^/(Tr. The numerical calculations are performed for the 
speckle potential described in section HII Bl and no fitting pa- 
rameter has been used. In addition, note that we have used 
a single realization of the random potential for each point in 
Fig. |4] Averaging over disorder turned out to have little im- 
portance, since the random potential is almost self-averaging 
in the range [— Ltf/4, +Ljp/4:], if CTr <C Ljp. 

Finally, we find from Eq. (fTSI l that the perturbative ap- 
proach that we have performed is valid whenever An <C no, 
i.e. whenever 



KVS(fTR/eo,0)« 



(29) 
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Note that this effect is more restrictive in the trap center where 
^0 is minimum. 



Figure 4: (color online) Amplitude of the fluctuations of the BEC 
density at the trap center, Arzc, in the combined harmonic plus ran- 
dom potential as a function of the ratio of the healing length to the 
correlation length of the disorder for several amplitudes of the ran- 
dom potential. The dots correspond to exact numerical results in the 
Gross-Pitaevskii approach [Eq. ©J and the lines show the theoretical 
prediction [Eq. (I28H. 



IV. CONCLUSION 

In summary, we have presented an analytical technique, 
based on the perturbation theory, to compute the static wave 
function of an interacting BEC subjected to a weak potential. 
This applies to the case where both the healing length of the 
BEC (^) and the spatial scale of the potential (ctr) are much 
smaller than the size of the system (L), but whithout restric- 
tion for the ratio ^/ctr. In particular, we have shown that when 
the healing length is larger than the space scale of the poten- 
tial, the BEC is sensitive to a smoothed potential which can 
be determined within our framework. 

Applying these results to the case of a ID random potential, 
we have shown that the wave function of a static interacting 
BEC is delocalized, similarly as in the TF regime 1 15]. This 
is confirmed by numerical calculations. The results of this 
analysis show that, for an interacting BEC at equilibrium, the 
larger the healing length, the smaller the perturbation induced 
by the disorder. It is worth noting that the conclusions of the 
present work hold for static BECs in the mean-field regime 
and when the interaction energy dominates over the kinetic 
energy in the absence of disorder, i.e. when the healing length 
is significantly smaller than the BEC half size (^ <C L). Going 
beyond the mean-field regime, it is interesting to study the 
interplay of interactions, disorder and kinetic energy in a Bose 
gas for interactions ranging from zero (where localization is 
expected) to the TF regime (where the BEC is delocalized as 
shown in this work). This question is addressed in Ref. l49ll . 

Finally, we note that the transport properties of a BEC can 
show significantly different physics. For instance, localiza- 
tion has been studied in matter-wave beams |33ll and in the 
expansion of an interacting BEC ifTsi [l6l llTll . In the lat- 
ter two cases, localization indeed does occur although non- 
negligible interactions can modify the usual picture of local- 
ization 
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